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Abstract: We calculate the Ohmic and Hall conductivities from linear response 
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1. Introduction 

With the proposal of the AdS/CFT conjecture [l| an outburst of activity in the 
relationship between gauge field theories and gravitational theories has appeared. 
The number of topics in which this gauge/gravity correspondence is shedding some 
light on is increasing significantly, one of the most important being the clarification 
of non-perturbative QCD-like theories 0. 

This regime of QCD is important to understand the phenomenology in the ex- 
periments currently carried at RHIC and shortly at LHC. Until the appearance of 
the AdS/CFT conjecture the only tools that could give some information about 
non-perturbative QCD were mainly the sum rules, chiral perturbation theory and 
lattice. This last lacks of a way to include real time dynamics and is not tractable 
analytically. 

The gravity dual of QCD is not known yet, but some advances have been made 
to approximate it. The better understood example (and also the original one) of the 
duality gives the description of M = 4 supersymmetric Yang-Mills ||, a theory that 
only has matter in the adjoint representation (which we call gluons). 

Karch and Katz found a way to include fundamental matter (quarks and their 
M = 2 superpartners) in the so-called quenched approximation [|J, in which the 
number of flavor (fundamental) degrees of freedom is much less than the number of 
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color (adjoint) degrees of freedom. This procedure does not consider backreaction of 
the quarks, thus corresponding to non-dynamical flavour. 

The work in [|, || considered systems where Nf Dg-branes were embedded in the 
background of N c ^> Nf Dp-branes in the presence of finite baryon density, modelled 
by a non-trivial profile for the temporal component of the £7(1) gauge field in the 
worldvolume of the probe Dg-branes. Dealing with such systems one frequently 
finds that the equations of motion are coupled under certain generic assumptions. 
For example, one can consider fluctuations moving with finite momentum and study 
the longitudinal component of £7(1) gauge field. At finite baryon density this channel 
is coupled with the scalar channel, which describes the fluctuations of the embedding 
profile of the probe Dg-branes in the Dp/Dq system J7|. If instead of a £7(1) baryonic 
gauge group one considers a ££7(2) isospin group one finds coupled equations in the 
transverse vector channel (fluctuations of the gauge field in the directions orthogonal 
to the fluctuation) ||. 

This coupling of the fields in the bulk of the gravitational theory translates via 
the AdS/CFT correspondence to a mixing of the dual operator on the boundary 
gauge field theory. A clear understanding of how to treat these mixing operators is 
essential in order to approach QCD phenomenology from a gravity dual, and only 
recently this subject has been studied in depth || [K]]. 

In the present paper we will deal with an analytic example with an already known 
solution in order to compare the proposed method of dealing with operator mixing. 
We will calculate the Ohmic and Hall conductivities in the presence of an external 
magnetic field for a plasma of gluons and quarks given by a Dp/Dq intersection. 
The solution was obtained in a macroscopic setup in [TTJ. The present calculation 
differs from that because we will exploit the Kubo formula. Thus we deal with linear 
response theory and will not need an external electric field. Similar calculations have 



been done in specific cases in [|12|, [13|, [Tj . 

This paper is organized as follows: in section we describe the gravitational 
setup we are going to use. In section |] we define the retarded Green's function from 
the on-shell boundary action for the case of mixed operators. In section ^ we calculate 
the Ohmic and Hall conductivities and finish in section |5| with some conclusions. 



2. Gravitational setup 

In this paper we will work in a general setup, but having in mind the embedding of 
Nf probe .Dg-branes in the background of a stack of N c ^> Nf coincident Dp-branes, 
where comparison with previous work is possible. In these systems there is a global 
U(Nf) symmetry on the .Dg-branes whose abelian center can be identified with a 
baryonic £7(1) symmetry, which we will exploit in order to study the system at finite 
baryon density. A recent review on the study of such systems can be found in 
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In this paper we will work with a metric of the form 

ds 2 = g Q0 (r)dxl + gn{r)d3% + g rr (r)dr 2 + g ee {r)<\Vt 2 n . (2.1) 

This metric is adapted to the pullback of a background of Dp-branes into the world- 
volume of Dq-branes in the quenched approximation. In that case the embedding 
profile of the probe Dg-branes in the transverse space is specified by a field ip(r) 
which enters in fl2.1p as 

g rr (r) = Gif(r) + ^(r) 2 G^\r, ^(r)), ; gee{r) = G$\r, 0(r)), 

with the components of the original 10- dimensional metric given by the stack of 
.Dp-branes, dQ^ is the metric of a unit n-sphere and we have minkowskian signature 
(i.e., #oo ( r ) < 0). Throughout this paper we will call the first three spatial coor- 
dinates 1 {xi, x 2 , x 3 , • • • } = {x, y, z, ■ ■ ■}. Our results are generalized easily to more 
general setups, for example the solution of a stack of Dp-br&nes without probes. 

We are dealing with finite baryon density so we are forced to assume the existence 
of a horizon at a radius Th 0- We will assume that the goo component of the 
metric vanishes linearly on the horizon as r — > r# (this can always be fulfilled by 
an appropriate redefinition of the radial variable). In order to have a well defined 
temperature one can easily show that the radial component has to diverge with the 
inverse of the distance to the horizon g rr {r ~ r#) ~ r 2 r rH w hh 7 r a constant. Our 
results are invariant under radial redefinitions and so the specific behaviour of the 
metric near the horizon is just a matter of convenience. 

In addition we will consider a background magnetic field, consequently we have 
a field strength tensor of the form 

Fab 

= {Sarho - daodbr) A' t (r) + {5 ay 5 bx - 5 ax 5 by ) B z . (2.2) 

The constant B z we will call a magnetic field, even when p > 3. 
Having in mind the DBI action, it is useful to define the matrix 

lab{r) = g a b{r) + F ab {r) (2.3) 

which is a function only of the radial variable appearing in the metric ( |2.1|) and has 
only diagonal componets apart from the two antisymmetric components coming from 
the electromagnetic tensor 70 = — 7 r o and j xy = — Jy x - 

1 notice that hereafter it is assumed that p > 2. 
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Notice also that the ju factors are all the same, but this is not the true for the 
inverse of the j a b matrix 



Too = -luf(r) ; Trr = G(r)f(r) ; f(r) = (r - r H )%(r 

y)0 _ Jrr y_L _ 7ii «s _ J_ 

TOOTrr + (TOr) 2 ' (Tm) 2 + (Try) 2 ' 7« ' 

rr _ TOO _ Or _. _ rO _ — TOr 

TOOTrr + (TOr) 2 ' TOOTrr + (TOr) 2 ' 



(t«) 2 + (t,,) 2 ' 



where 7^, 7^, 7or, 7»y, £(r) and G(r) are regular at the horizon and "y ab j ca = & c - 7~ L ~ L 
stands for 7 XX = j yy and j zz for all the other p — 2 flat-spatial diagonal components 
of the 7 matrix. 

2.1 Background field solution 

To write a solution for the background gauge field components we will be focusing 
on the DBI action 

Sdbi = -N f T Dq J d p x dr dfi n e^y^, (2.4) 

where 7 = det( 7afc ) = gfr 2 g% g (g 00 g rr + (2na') 2 A' t 2 ) ((g u ) 2 + (27ta') 2 B 2 ) and dfi n is 
the volume form of a unit n-sphere. Notice that we are not considering Wess Zumino 
terms, which are not needed in the present case. Its presence is discussed further in 
FT, [TJ], [RJ . The equation of motion one reads from ( |2.4| ) for the background gauge 



field is 



SSDBI _ AT rr, fn /\-d, /—nr n Q 



^ - N f T Dq (27ra')e-^l° r = (2.5) 

with Q n the volume of the unit n-sphere. As all the angular dependence is encoded 
in the volume form, we can identify the constant n q with the density of quarks in 
the field theory || . We can use this constant of motion n q to find a solution for the 
gauge field component 



A' t (r) = n q J = , (2.6) 

9 V APe-^r 2 (ti + (2na>) 2 B 2 ) <fc + (2na') 2 n 2 q ' 

with Af = NfT Dq Q n (27va') 2 . The holographic dictionary tells us that the value of 
this field at the boundary is related to the chemical potential, /i, for the quarks, 
whereas at the horizon one has to impose vanishing of the gauge field component to 
avoid having a singular one- form ||, thus we find for the chemical potential 



/ / ~ 9 00 9 rr j 

"' xl ^ 2 e- 2 ^- 2 (g 2 + (2 7 ra') 2 B 2 )g^ + (2na') 2 n 2 
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For the embedding profile of the probe D7-branes ip(r), the equation of motion 
one finds is a highly non-linear one. The solutions depend non-analytically on two 
parameters ipQ (the angle at which the branes enter the black hole) and n q |J. 

3. Holographic retarded correlator for mixed operators 

Mixed operators arise naturally in holographic models in the form of coupled dif- 
ferential equations for the fluctuations of the background fields. On occasions these 
couplings involve non-physical fields, as is the case of the temporal and longitudinal 
fluctuations of a one-form, which can be combined in terms of a gauge invariant 
quantity to decouple the system [[17]]. On the other hand, at times the coupled equa- 
tions appearing already involve physical gauge invariant fields. When this occurs the 
operators get mixed in the IR of the gauge theory. This mixing can be propagated 
to the UV. 

There are several examples in the literature where one finds this mixing of op- 
erators. In particular it appears in the study of the D3/D7 system at finite baryon 
density as a coupling between the longitudinal and scalar perturbations 0. A sim- 
ilar system, but with finite isospin density instead of finite baryon density, also has 
such a coupling in the transverse channel, giving rise to a superfluid behaviour of 
the system ||. These issues present some subtleties on the study of these systems, 
for example, if one wants to study the quasinormal modes of a system with operator 
mixing one has be particularly careful, as shown for example in M, [HJ. 

We will bring here a simplified discussion of the holographic prescription to study 
this particular kind of systems. The linearized equations of motion we want to solve 
can be derived from the action 2 

S = ~y y"d p xdre-*v^7^<l T -r • 0„$ , (3.1) 

where T is a matrix which couples the different fields (given by the column vector 
x p , r)) in our bulk theory. The form of this matrix is responsible for the mixing 

of operators at the boundary. 

Assuming the 7 components are of the form (|2.3j ), the Euler-Lagrange equations 

after expressing the fields in the mode representation x p , r) = e~ ti < U)t ~ qXpS, Q(r) 

read 

d r {e~*^y r T s ■ &)+[icud r (e-^v^T "^) - e"*V=7 (u; 2 7 00 + gV) T s ] •$ = 0, 

(3.2) 

with TsfA) the (anti) symmetric part 3 of T. 

2 note that this action for the fluctuations is just a Maxwell-like one, so this discussion is inde- 
pendent of the fact that it originally comes from a DBI action. 

3 we follow the convention Tg = T+ 2 ' r and similarly for the antysymmetric part. 
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One can express the on-shell boundary action from Q3.1|) as 



S 



bou 



K 

Y 



d p x e 



- 7 $ J 



Y r r s ■ <9 r $ - 7 ur r A • a $ 



„ 1 A 



(3.3) 



in 



where we evaluate it at a cutoff r = A. Reexpressing it as the product of a profile and 
a source x p ,r) = e - l (ut-qx P ) ] F k (r) ■ ip(k) with F k (A) = I; the Fourier transformed 
action reads 



Si 



bou 



K 

T 



d p k _ A 

(p T (-k) ■ T k ■ ip(k) 



(2tt 



(3.4) 



The gauge/gravity conjecture relates this on-shell boundary action with the re- 
tarded Green's function by 



G R (k) = - lim ^(A), 



A^oo 



which for the case of unmixed operators (i.e., a diagonal T) coincides with the 
prescription given in |18|], which has been derived also via the Schwinger-Keldish 
formalism [19, ^Of. One can show that G R (—k) = G R (k)* = G A (k), see appendix A 
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When dealing with unmixed operators it is known that one can define a conserved 

flux 

d r lmG R (k) = 0. 
This is promoted to the condition 

d r (T k - H) = F? k ■ d r (e~t^lY r Ts ■ F' k ) - d r [e^^f r (F^)' ■ T s ) ■ F k = , 

(3.5) 

where we used T k = T- k and the last equality is guaranteed by the equations of 
motion. This conserved flux can also be seen as the conservation of a Noether 



current due to a U(l) symmetry of the bilinear action [10 



4. Fluctuations and the conductivity matrix 

We will now introduce perturbations to our system. Given the existence of a non 
vanishing background field B z the most general perturbation we can write is ex- 
pressed as a mode e - l \ut-<imx-<izz) _ Having the Hall effect in mind we want to make 
perturbations of the electric field in the directions perpendicular to the background 
magnetic field. For this purpose we will restrict to the simpler case of a fluctuation 
parallel to the background magnetic field (q x = 0), in which the possible perturba- 
tions appearing can be classified according to their transformation under the little 
group SO(p-l) 0. 
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The Kubo formula we will use to calculate the conductivity matrix requires the 
evaluation of the transverse channel given by the perturbations Ap^ z (t, z, r). The 
presence of a magnetic field given by a gauge potential A x (y) couples the fields 
Ax and A y , leaving the other transverse components independent. Including the 
perturbations, we shall deal with a gauge field of the following form 



A„ 



(A t (r), B z y + A x (t, z, r),A y (t, z, r), 0, • • ■ , 0) 



(4.1) 



which amounts to doing perturbations on top of a macroscopic magnetic field. One 
can find the solution for the other transverse perturbations Ap^{ X) y }Z } from these two 
by taking the B z — > limit. 

To study the linearized equations of motion one should expand the DBI action 
to second order in fluctuations. One can do this and then compare to the general 
expression (|3T ), finding 

T 



(4.2) 



These equations can be solved in the hydrodynamic limit 4 , as shown in the 
appendix. The final solution for $ in this limit is 



$0) 



I 



% UJ 



C 



l UJ 



rfUrf^ /0 -1 

_ 1 



c 2 l ±X l rr 



•<p(k) 
(4.3) 



= F k (r)-tp(k), 

where C*2 is chosen such that Fk(A) = I, (p(k) = $(A) and C is the constant matrix 



C 



"7 



(4.4) 

is conserved at the given 



It can be shown that the antihermitian part of the flux 
order in u> and q z . 

Following the formalism described in section [3] we can read off the retarded 
transverse Green's function as 



Gf(k) = - lim F k (A) = jV lim e~*^ \T S ■ F' k + iuJ^ r T A ] 

A^oo — "> 



-I UJ 



AfC H , 



where the dots express higher order in the frequency and momentrum. From the 
former equation we can evaluate the matrix of conductivities as 



o\uj 



■G*(uj=\q z \) 



UJ 



jVC + 0(uj). 



(4.5) 



4 in this limit the equations decouple in the appropriate basis, but this fact does not affect our 
discussion on operator mixing. 
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As we are working in the hydrodynamic regime we can only study the DC con- 
ductivity cr(0), given by 



(7(0) 



a 



■>:y 



- <T X y <T XX 



where 



a, 



Me 



' \An /00 7 rr 7" L± 



cr 



'7/ 



'7/ 



(4.6) 



(4.7) 



This result depends only on the bilinear action ( |3.1| ) with coefficients ( |4.2| ) and 
the presence of a suitable black hole, so this formula can apply with minor modi- 
fications to other setups apart from Dp/Dq intersections. Due to the linearization 
of the perturbations the importance of the original DBI action in these results gets 
reduced. 

Additional transverse components (appearing when p > 3) will give a diagonal 
conductivity given by the B z — > limit of fl4.7|) , i.e. 



a 



Me 



7(s z 



,o)7°V r 7" 



This is the result obtained in [|Tj| [2T| . 

When the dependence on the charge density is explicitly needed one has to plug 
a solution for the background gauge field A t . In this case the form of the original 
action becomes relevant. For example, plugging the results (|2.5|) and ( |2.6| ) in (|4.7|) 
one obtains for the DC conductivity and the Hall coefficient as originally found in 

H 



AT- 



Or 



+ 



7l + (2<Ka>yBl ( 7 2 + [2^a'fBlY 



a 



(2na 



A2 



n q B z 



in 



hii(r H )) 2 + (2tt«0 2 S z 2 
(4.8) 

The insertion of a background electric field in our setup would have complicated 
the analysis performed. An external electric field induces the presence of a shell 
outside the horizon (described by a radius related to the modulus of the electric 
field) on which all the fields must be regularized. In a black hole setup this induces a 
current on the boundary that can be analyzed to get the values for the conductivity 
and the Hall coefficient, whose zero electric field limit is given by ( |4.8j ). 

As noted in JIT], [14], IB], adding a 9 term to the original action has as an effect 
on the present result a shift proportional to 9 in the charge density and in the Hall 
coefficient. 



5. Conclusions 

We used a generalized prescription to obtain the Green's function in the case of 
operators mixing on the boundary to calculate the conductivity and the Hall coef- 
ficient of a plasma with a gravity dual given by a system of Nf -C N c Dg-branes 
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in the background of N c Dp-branes. The solution to this problem in a macroscopic 



setup was found in [ 11 , g2| , where the matrix of conductivities was computed via the 



relation between induced currents on the boundary and external electric fields 

The present calculation complements the previous, macroscopic approach with the 
microscopic procedure, making use of the Kubo formula ( f4.5|) . Agreement between 
both procedures is found. 

Both methods are based in the regularisation of the fields at a given radius, set 
by the external electric field. When this electric field is set to zero the radius at 
which one has to regularise coincides with the radius of the horizon. The transport 
phenomena information one can extract from this kind of systems is encoded uniquely 
by this surface, and it appears naturally by imposing regular behaviour. 

An important point given in |TT| is that the expressions (|4.8j), evaluated at the 



radius r* instead of the horizon, would give the result for the conductivity and the 
Hall coefficient in the presence of an electric field E creating a singular shell at 
precisely that radius r*(E). We suggest that this still holds in the more general 



result (|4.7| ). In |23| the conductivity was calculated in a setup with just an external 



electric field and this suggestion was found to be true. 
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A. Solving the equations of motion 

The equations of motion (|3.2| ) can be written as 



V + log' [ c - V^-T] <&' - 7 ^ * 

xy Or 

+u\og' [ e ^v^7 x V r ] -^ry ■ * = 0, (A.l) 

with a 2 the second Pauli matrix. This can be solved by shifting to the basis of 
polarizations 
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in which the equation becomes 
+ log' 



i re~^r 



xy Or 

-a; log' [^V^T*] ^j^a 3 • M>, 



(A.2) 



where a 3 is the third (and thus diagonal) Pauli matrix. Thus, these are two uncoupled 
equations for the components of 

The horizon (r = r#) is a regular singular point, therefore we can perform the 
usual Frobenius study to find its singular behaviour. The indices at the horizon are 



4^f 



and we will consider just the negative sign to regularize the function in that point 

A± f 7 * •A-± ) reg ■ 

A.l Hydrodynamic solution 

We can make an expansion at low to and q. To do this we set 

(u, q z ) -»• Xhyd(uJ, q z ) , 

A ± - P- [AS + \ h y d A { l ] ) = A { ° ] + \ hyd (A { 1 ] + r/_^ 0) log/) . 

One finds that A^± is a constant and an equation for A*± 



which can be solved as follows 
A ( l\r) - 



c 2 



(A.3) 



(A.4) 



where the constant C2 is the same for both solutions. 
We can now express the solution for A x>y as 



|(0) 



2 



+ ue -<t>. f^^xy^r 



1T"T' -^— 2 — 



c 2 



+ 0(\L d ), (A.5) 



Ay — 



Af - A 

2i 



(o) 



+ 



c 2 



j(0) + (1 («) 
— xy 0r*+ 

+ 0(X 2 hyd ) . (A.6) 



In order to fix the integration constants A±* we normalize the fields A x , y to the 
value at the cut-off r = A. For this it suffices to take A*± = A% ± iA y and C 2 such 
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that the second term in the solution vanishes at the cut-off, and from this we read 

C x + iwe-^v/^rf f 7 0r ^ 



Jc 2 



r C y - iue-* y /=i'f yr ) 0r d£ 



+ o(\ 2 



hydJ •> 



g — 4> ^J^—ry ry 1 1 ry T V 



(A.7) 
(A.8) 



The regularizing constants C x and C y are fixed by imposing regularity of the 
function Ax} y ,reg{r) = Ax) y {r) — log fA^y on the horizon. For A x we obtain after 
derivation the condition 



C x + lue-^^i^^A^ iuj f 



finite , 



4?rT / 

and we finally obtain (using the fact that e - ^^/— 77^7°'* is finite on the horizon) 



C r = —ioJ e 



C y = —iu e 



"7 



7°V r 7^ + l xy l° r A^ 



y/-'f r r l ii ^ - l xy l° r A^ 



r—>rn 



(A.9) 
(A.IO) 
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